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ABSTRACT 

a7 

A  code  has  been  written  to  use  the  algebraic  computer 
system  MACSYMA  to  generate  systematically  the  in¬ 
finitesimal  similarity  groups  corresponding  to  systems 
of  quasi-linear  partial  differential  equations.  The 
infinitesimal  similarity  groups  can  be  used  to  find 
exact  solutions  of  the  partial  differential  equations. 
In  an  example  from  fluid  mechanics  the  similarity 


method  using  the  computer  code  reproduces  immediately  a 
solution  obtained  from  dimensional  analysis. 


1.  Introduction 

In  recent  years  Lie  infinitesimal  contact  transformation 
group  theory  has  become  a  widely  used  tool  for  systematic  ex¬ 
plorations  of  similarity  solutions  of  partial  differential 
equations  [1] .  Though  the  implementation  of  this  method  in¬ 
volves  mostly  standard  manipulations  (algebra  and  chain  deri¬ 
vatives)  ,  in  practice  these  become  extremely  involved  for  all 
but  the  simplest  of  equations. 

Here  we  introduce  the  use  of  the  algebraic  computing 
system  MACSYMA  to  facilitate  these  calculations.  Specifically, 
MACSYMA  is  used  to  calculate  systematically  the  generators  of 
the  infinitesimal  group  under  which  the  considered  equations 
are  invariant.  As  far  as  the  operational  side  of  the  similarity 
method  is  concerned,  this  is  the  most  involved  and  tedious  step. 
Once  the  general  similarity  group  is  found,  a  (hopefully  non¬ 
trivial)  subgroup  is  found  which  leaves  invariant  boundary 
curves  and  boundary  conditions.  Finally,  once  this  subgroup 
is  found,  its  invariants  and  consequently  the  similarity  forms 
of  the  solutions  of  the  partial  differential  equations  can  be 
found. 

A  primary  advantage  of  the  similarity  method  is  that  it 
is  one  of  the  few  general  techniques  for  obtaining  exact  (non¬ 
linear)  solutions  of  partial  differential  equations.  A  primary 
disadvantage  of  the  similarity  method  is  that  the  solution  found 
may  satisfy  only  a  very  restricted  set  of  initial  and  boundary 
conditions . 
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In  general,  the  existence  of  a  one  parameter  group  of 
transformations  leads  to  a  reduction  by  one  in  the  number  of 
independent  variables.  Thus,  in  the  most  encountered  case  of 
two  independent  variables,  the  reduction  leaves  us  with  an 
ordinary  differential  equation. 

We  remark  here  that  the  similarity  method  contains  as  a 
special  case  those  solutions  which  can  be  obtained  via  dimen¬ 
sional  analysis.  That  is,  if  an  equation  has  a  solution  which 
can  be  obtained  from  dimensional  analysis,  then  the  present 
method  will  reproduce  immediately  that  special  similarity 
solution,  and  possibly  more  general  similarity  solutions  as 
well. 


2 .  Sketch  of  the  Similarity  Method 

In  this  section  we  present  a  brief  summary  of  the  similar¬ 
ity  method.  The  similarity  transformation  solver  (STS)  uses 
MACSYMA  to  construct  the  similarity  transformations  which  are 
admitted  by  a  set  of  partial  differential  equations  of  the 
form  (sums  over  j,k,  and  1) 


aij 


at 


u.  +  b.  .  rr-  u.  +  c 

-i  i  t  at  t 


2  "j 


..  u 


+  d.  , 


ijkl  ax^Xj^  uj 


ei  ' 


where  u^  «  Uj (x,t)  are  solutions,  and  a,b,c,d,  and  e  are 
arbitrary,  specified  functions  of  (x,t,u)  . 


The  one  parameter 


infinitesimal  transformations  considered  are 

x'  =  x  +  eX(x,t,u)  (2a) 

t'  =  t  +  eT(x,t,u)  (2b) 

u'  =  u  +  eU(x,t,u)  .  (2c) 

MACSYMA  is  used  to  determine  the  differential  equations, 
called  the  determining  equations,  that  are  satisfied  by 
the  infinitesimal  groups  j£,T,ij.  Technically,  this  is  achieved 
by  requiring  that  the  partial  differential  equations  (1)  be 

2  4-  4 

invariant  to  order  £  in  terms  of  the  new  variables  (x',t',u') 
of  Eqs.  (2).  Thus  the  annuling  of  0(e)  terms  becomes  an  in¬ 
variance  condition. 

Let  M[u]  =0  be  a  system  of  partial  differential 
equations  in  the  old  variables,  and  let  M'[u']  =  0  be  the 
same  system  in  the  new  variables,  where  is  obtained 

from  M[u]  by  replacing  everywhere  (x,t,u)  ■+•  (x',t',u') 

The  application  of  the  transformations  (2)  to  M'fu'J  shows 
that  M'[u']  is  of  the  form  M'[u']  =  M[u]  +  eH  +  0(e2), 
where  H  «  &(x,t,u,X,T,i5)  .  We  make  the  system  invariant  to 
the  infinitesimal  transformations  (2)  by  requiring  that  H  =  0. 
Suppose  we  denote  by  {Sm)  »  for  an  appropriate  set  of  values  m, 
all  possible  first  or  second  derivatives  of  u^  with  respect 
to  x^  and/or  t  .  Then  the  expression  $  is  a  sum  of  terms 
involving  the  Sm's,  with  coefficients  functions  of  (X,T,6) 
and  their  derivatives  with  respect  to  (x.  ,t,u. )  .  After  the 

K  j 

original  system  of  partial  differential  equations  is  used  to 
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'i*  w.  ■pi-i,.  ,  i -  Juwmtwwwinp 


guarantee  that  all  the  Sm's  that  appear  in  H  are  independent, 

H  is  made  to  vanish  by  requiring  that  all  the  coefficients  of 

the  Sm's  vanish.  This,  as  explained  in  [1],  is  a  sufficient 

but  not  a  necessary  condition.  The  coefficients  of  the  S  's 

m 

are  differential  equations  satisfied  by  the  infinitesimal 
generators  (X,T,U);  these  are  the  determining  equations  mentioned 
earlier.  In  practice  the  determining  equations  usually  can  be 
solved  easily,  as  they  are  an  overdetermined  set  of  linear 
equations.  A  knowledge  of  the  infinitesimal  transformations 
enables  one  to  construct  a  continuous  transformation  or  in  turn 
to  find  the  invariants  of  the  group  that  serve  as  the  new  de¬ 
pendent  and  independent  variable,  it  is  in  these  new  variables 
that  Eqs.  (1)  exhibit  their  symmetry  and  result  in  a  re¬ 
duction  by  one  in  the  number  of  independent  variables. 


To  illustrate  the  desirability  of  using  MACSYMA  to  calcu¬ 


late  the  determining  equations,  suppose  a  partial  differential 

equation  contains  a  derivative  (assume  scalar  u  and  x) . 

In  terms  of  the  new  variables,  the  transformed  partial  differential 

3u' 


equation  contains  a  term 


dX 


But  from  Eq.  (2c) , 


Ah  *  o 

-  =  -  [u (x , t)  +  eU (x, t,u) ] 

3x'  3x' 

_  3u  3x  .  3u  3t  .  I  3U  3  x  3U  3t 


+ 


3U 

3u 


3u  3x 
**  9x> 


(3) 
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1 


Now  for 


3x 

3x' 


and 


3 1 
3x' 


use  Eqs.  (2a-b)  to  obtain 


3x  = 
3x' 


1  - 


e 


and 


(4) 


(5) 


By  substituting  Eqs.  (4)  and  (5)  into  Eq.  (3)  and 
simplifying,  we  obtain 


3u'  3u  [30  f3U  3X1  3u  3T  3u  3X  [3ul 2 

3x,  3x  [_3x  [3u  3xJ  3x  3x  3t  3u  [3xJ 

-HI? It)  +°<e2>  •  '6> 

32u 1 

Correspondingly,  any  second  derivatives  such  as  - become 

3x' 

very  involved.  The  system  of  partial  differential  equations 
(1)  is  written  in  a  completely  vectorized  form.  Since 
MACSYMA  knows  the  chain  rule  and  can  simplify  and  factor,  it 
is  quite  feasible  for  MACSYMA  to  reliably  transform  compli¬ 
cated  systems  of  partial  differential  equations  in  several  in¬ 
dependent  variables. 


3 .  Programming  Considerations 

In  this  section  we  describe  the  implemention  of  the  STS 
on  the  MACSYMA  system.  The  basic  mode  of  operation  is  that  a 
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reference  version  of  the  STS  is  modified  as  required  with  the 
TECO  editor  to  solve  the  problem  at  hand,  and  then  the  modified 
version  is  executed  while  in  MACSYM  using  the  BATCH  command. 

Once  program  execution  is  initiated  with  the  BATCH  command, 
no  further  user  interaction  is  required. 

Input  procedures  for  the  STS  are  described  as  follows. 

Once  a  user  obtains  a  reference  copy  of  the  STS  program  in  his 
user  space,  a  listing  of  the  program  should  be  made  while  in 
TECO.  In  the  program  listing  the  portions  of  the  code  which 
must  be  changed  from  one  problem  to  another  are  preceded  by  a 
single  asterisk  and  a  blank  line  and  are  followed  by  a  blank 
line  and  a  double  asterisk.  There  are  three  places  in  the 
STS  program  where  the  user  must  modify  the  text  to  solve  a 
new  problem.  First,  a  title  should  be  specified  which  identifies 
the  system  of  partial  differential  equations  being  considered. 
Second,  the  user  must  specify  the  number  of  equations  (i.e., 
the  number  of  dependent  variables  u^) ,  NEQ,  the  spatial 
dimension  of  the  problem,  NX  (NX  1) ,  and  whether  or  not 
there  are  time  derivatives  (NT  =  0,  no  time  derivatives; 

NT  =  1,  there  are  time  derivatives) .  Finally,  the  user  must 


specify  the  nonzero  coefficients  a..,  b..,  c...,  d.  and 

13  lj  ljjc  ljjvj. 

e^  or  Eq.  (1) .  Before  assigning  the  subscript  i  to  the 
coefficients,  the  equations  should  be  ranked  in  ascending 


order  according  to  the  number  of  terms  in  the  equations.  For 


example,  the  equation  with  the  fewest  terms  goes  first,  etc. 


After  these  changes  are  made  in  the  reference  STS  program,  the 


-7- 


updated  version  is  stored  with  TECO  as  a  new  file,  to  be 
excuted  later  with  BATCH.  In  this  way  a  user  might  build  a 
library  of  STS  programs  appropriate  to  different  systems 
of  equations. 

The  output  of  the  STS  is  in  two  stages.  In  the  first 
stage  the  BATCH  command  causes  listing  and  execution  of  the 
code.  During  this  stage  some  of  the  intermediate  calculations 
are  printed  out,  so  that  the  progress  of  the  calculation  can 
be  inspected  as  it  proceeds.  This  portion  of  the  computing 
time  takes  about  10  minutes ,  nearly  independent  of  the  system 
of  equations  considered.  The  first  intermediate  printout  is 
a  listing  of  the  functional  dependencies  of  the  functions  in¬ 
volved  in  the  infinitesimal  transformations  (2) .  The  STS 

program  is  vectorized,  so  the  variables  (x,  , t ,u . ,X.  ,T,U . )  are 

Jc  1  k  j 

replaced  with  the  program  variables  (X [K] ,T,U [ j] ,FX [K] ,FT, 

FU [ J ] ) ,  respectively.  The  second  intermediate  printout  is  a 
listing  of  the  system  of  equations  to  be  solved.  This  serves 
as  a  check  that  the  equations  were  ordered  properly  and  that 
the  coefficients  were  entered  correctly.  Next  the  linearized- 
in-epsilon  versions  of  the  new  coefficients  of  the  transformed 
partial  differential  equations  are  displayed.  Recall  that  the 
original  system  of  equations  must  be  used  to  guarantee  that 
the  S^'s  that  appear  in  H  are  independent.  The  fourth 
intermediate  printout  lists  NEQ  of  the  Sm's  which  the 
STS  has  chosen  to  eliminate  from  H  by  using  the  original 
system  of  equations.  The  solutions  of  the  original  equations 


for  the  NEQ  S  's  is  then  listed  in  the  last  intermediate 
m 

printout . 

The  second  stage  of  output  in  the  STS  is  the  displaying 
of  the  components  of  H,  from  which  the  determining  equations 
are  derived.  Each  component  of  H  is  factored  conveniently 
with  respect  to  all  the  S^'s;  t‘he  determining  equations  can 
be  read  off  immediately  as  the  coefficients  of  the  independent 
terms  involving  the  Sm's*  The  amount  of  time  spent  computing 
H  depends  significantly  on  the  system  of  equations  solved,  as 
well  as  on  the  current  load  on  the  MACSYMA  system.  A  very 
rough  estimate  of  the  computing  time  to  calculate  2  is 

T  (min)  =  NEQ  (NX+1)  (NX/2  +  2).  (7) 


4 .  An  Example 

As  an  example  of  the  use  of  the  STS,  consider  the  classical 
Blasius  problem  from  fluid  mechanics;  solve  the  problem  of  two- 
dimensional,  steady  state,  incompressible,  viscous,  laminar  flow 
in  a  boundary  layer.  The  equations  to  be  solved  are 


5y  vy 


g  g  9  2V 

vxHVx  +  vy  5y  vx  =  v  JJT 


(8a) 

(8b) 


where  v  =  (v  ,v  ,0)  ic  the  fluid  velocity,  and  v  is  the 
x  y 

coefficient  of  viscosity.  A  semi-infinite  plate  exists  in  the 
x  -  z  plane  for  x  >_  0  ,  and  far  from  the  plate,  y-*-±®  ,  vx 
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assumes  a  constant  value. 


The  viscous  boundary 


conditions  are: 


u 


o 


vx(y  =  ±°°)  =  uQ  (9a) 

v„(y  =  0)  =  v  (y  =  0)  =  0  .  (9b) 

x  y 

Using  arguments  of  dimensional  analysis,  one  concludes  that 
(Landau  and  Lifshitz  [2])  the  new  independent  similarity 
variable  is 


c  =  y  /uQ/x 

(10a) 

and  that  v  and  v  are  of  the  forms 
x  y 

-  °o  fl(c> 

(10b) 

vy  -  (%]  f2  (t>  ' 

(10c) 

Let  us  now  obtain  these  same  results  as  well  as  the  gen¬ 
eral  infinitesimal  similarity  group  by  means  of  the  STS.  In 
a  program  listing  of  the  STS  we  first  insert  in  the  appropriate 
place  a  title  for  problem  being  solved.  Next  we  insert  on 
separate  lines  in  the  STS: 


Finally  we  insert  the  nonzero  coefficients: 


C [1, 1/ 1] : 1$ 

011,2,21:1$ 

C  [2, 1, 1]  :U  [1]  $  (12) 

C[2,l,2] :U [2] $ 

D [2, 1,2, 2] :  -  NU$  , 


where  NU  is  the  coefficient  of  viscosity,  v  .  After  this 
version  of  the  STS  is  stored  in  TECO  and  executed  in  MACSYM 
with  BATCH,  the  printouts  begin. 

The  components  of  H  are  printed  out,  and  the  determining 
equations  can  be  read-off  by  inspection  (see  the  Appendix) . 

They  are  (letting  U[l]  u,  U[2]  -*•  v,  FU[1]  -*■  U,  and  FU[2]  -*•  V) 


X  =  X  =  X  =  X  =  Y  =  Y 
v  u  y  w  w  v 


=  U  =  u  =  Y  =  0 
w  v  u 


V  -  Y  =  V  +  U  =  0 
u  x  y  x 


(13) 


Y  -  V  -  X  +  U  =  U  -  uX  +  2uY  =  0 
y  v  x  u  x  y 


V  <“Uuu  +  2Yuy)  +  2vYu  *  -  vUyy  +  vUy  +  uUx  =  0 


v  (Y  -2U  )  +  vY  -  uY  +  V  =  0  . 

yy  uy  y  x 


This  set  of  differential  equations  for  X,Y,U,  and  V  as 
functions  of  x,y,u  and  v  is  overdetermined. 
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The  integrations  involved  is  solving  Eqs.  (13)  are  straight¬ 
forward.  The  result  is  that  the  general  infinitesimal  similarity 


group  corresponding  to  Eqs.  (8)  is 

X(x,y,u,v)  =  (a  +  2b)x  +  c 
Y(x,y,u,v)  =  by  +  g  (x) 

U (x,y ,u,v)  =  au 
V(x,y,u,v)  =  ugx  -  bv 

where  a,b,  and  c  are  arbitrary  constants,  and  g(x)  is 
an  arbitrary  function  of  x  alone. 

Now  we  find  the  subgroup  of  (12)  that  leaves  invariant 
the  boundary  curves  (9) : 

Invariance  of  y  =  0=*y'  =  0  =»  Y  (x,y  =  0,u,v)  =  0 
**  g(x)  =0  . 

From  the  transformations  (2a) ,  the  boundary  curves  y  ■  ±  «*> 
are  trivially  invariant. 

Invariance  of  x  =  0*»x'  =  0*»X  (x=  0,y,u,v)  =  0 
♦  c  *  0 

Invariance  of  x  >  0  ■*  x'  >  0  *  X  (x  >  0,y,u,v)  >  0 
2b  +  a  >  0  . 

Equations  (15)  will  restrict  the  group  (14)  so  that  the 
boundary  curves  are  left  invariant.  Now  we  further  restrict 


(14) 


(15a) 


(15b) 


(15c) 
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the  group  (14)  so  that  the  boundary  conditions  are  invariant: 
u (y  =  0)  =  0  ■*  u'  (y  =  0)  =  0  U (x,y  =  0,  u=0,v)  =  0 

♦  no  restriction  on  a  .  (16a) 

Similarly,  v(y  =  0)  =  0  is  trivially  satisfied 

no  restriction  on  b  .  (16b) 

u(y  =  ±  ®)  =  uQ  ■*  U  (x,y  =  ±  »,  u  =  uq,v)  =  0 

*  a  =  0  (16c) 

By  combining  Eqs.  (15-16)  the  final  subgroup  of  Eqs. 

(7-8)  that  leaves  invariant  the  boundary  data  is 

X(x,y,u,v)  =  2bx 

Y(x,y,u,v)  =  by  (17) 

U(x,y,u,v)  *  0 
V (x,y , u,v)  *  -  bv. 

The  first  invariant  surface  condition  is 


dx  _  dy  _  du 
X  Y  u  ' 


or. 


dx  dy  _  du 

2bx  by  o 


Integration  of  the  first  equality  in  (18)  is 


C  *  ^ /  /x  *  constant. 


Integration  of  the  second  equality  in  (18) 


(18) 


(19a) 


u  »  F1(C) 
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(19b) 


The  second  invariant  surface  condition  is 


dx  _  dy  _  dv 
X  Y  ”  V 


or. 


dx 

2fc>x 


dy  dv 
Ey  E v  * 


Integration  of  the  second  equality  in  (20)  gives 


(20) 


v  «  f2(c)  /y  •  <19c> 

Equations  (19)  are  the  same  results  as  Eqs.  (10)  obtained 
by  Landau  and  Lifshitz.  Furthermore,  the  infinitesimal  group 
(14)  is  the  general  infinitesimal  group  corresponding  to  Eqs. 

(8) .  Thus  the  similarity  solution  corresponding  to  any  other 
boundary  conditions  can  be  found  by  analyzing  the  group  (14) . 
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Appendix 


In  this  appendix  we  illustrate  how  the  determining 
equations  are  obtained  from  a  printout  of  the  functions  H  . 

As  an  example  we  consider  the  system  of  partial  differential 
equations  (8),  so  there  are  two  components  H^i  =  1  and  2. 

The  STS  factors  each  IL  to  produce  a  polynomial  in  the 
Sm's,  where  the  factoring  is  performed  first  with  respect  to 
Sjy  then  with  respect  to  S2,  and  so  on.  That  is,  HL  first 
is  factored  to  produce  a  polynomial  in  S^  .  Then  the  coef¬ 
ficient  of  each  power  of  S^  is  factored  to  produce  a  poly¬ 
nomial  in  S2,  etc.  In  this  way,  once  a  particular  product 
of  powers  of  the  Sm's  appears  in  ,  that  same  combination 
of  the  Sm's  never  appears  again.  Consequently  the  coef¬ 
ficients  of  the  independent  terms  of  HL  (i.e.,  the  deter¬ 
mining  equations)  can  be  identified  by  inspection. 

The  printout  of  H  corresponding  to  the  partial  differential 


COMPONENT  NUMBER  1  OF  THE  EXPRESSION  H  PROM  WHICH  THE 
DETERMINING  EQUATIONS  ARE  DETERMINED  IS: 


d  d  2  d  d  d  d 

ZEROP  =•  << - PX  - FX  >*S  +  < - FX  - FU  - FX  + - FU  >*S 

dU  2  dU  14  dX  2  dU  2  dX  1  dU  1 
2  1  2  2  11 


d  d  d  d 

+  (( - FX  FX  )*S  - - FX  + FU  >*S 

dU  2  dU  1  13  dX  2  dU  2  7 

2  1  11 

d  d  d  d 

+  < - FX  + - FU  >*S  + - FU  + - FU  >*EPSLON 

dX  1  dU  1  13  dX  2  dX  1 

2  2  2  1 


COMPONENT  NUMBER  2  OF  THE  EXPRESSION  H  FROM  WHICH  THE 
DETERMINING  EQUATIONS  ARE  DETERMINED  IS: 


d  d  d 

ZEROP  =*  <<-  2*< -  FX  >«NU«S  +  2*< -  FX  >*NU*S  +  2*< -  FX  >*NU>*S 

dU  1  4  dU  1  7  dX  1 

2  12 


2  2  2 
d  3d  d 

+  < - FX  )*NU*S  +  << - FX  -  2*  < - FX  >>*NU*S 

21  4  22  dUdUl  7 

dU  dU  12 

2  2 

2  2 

d  d  d  d  2 

+  <-  2* < - FX  ) - FU  >*NU  -  2*U  *( - FX  )  +  U  #< - FX  >>*S 

dU  dX  1  2  1  1  dU  2  2  dU  1  4 

22  dU  2  2 

2 

2  2 
d  d  2 

+  <<-  2*< - FX  )  + - FX  )*NU*S 

dU  dU  2  2  1  7 

12  dU 

1 


18 


+  <<-  2*< 


2  2  2 
d  d  d 

- FX  )  +  2*< - FX  )  +  2*< - FU  )  >*NU  -  l 

dU  dX  2  dU  dx  1  dU  dU  1 

2  2  12  12 

d  d  d 

+  2*U  *< -  FX  >>*S  -  U  *< -  FX  >*S  +  ( -  FX  >*NU*S 

1  dU  2  7  1  dU  1  13  dU  1  14 

1  2  2 

2  2 

d  d  d  d 

+  < - FX  +  2* < - FU  )  )*NU  +  2*U  *< - FX  >  -  U  *< - FX  ) 

2  1  dU  dX  1  1  dX  2  2  dX  1 

dX  2  2  2  2 

2 

2 

d  d  d  3 

-  U  *< -  FX  )  -  U  *( -  FU  )  +  FU  >*S  +  < -  FX  >*NU*S 

ldXl  2  dU  1  14  22  7 

12  dU 

1 

2  2 

d  d  d  2 

+  <<2*< - FX  ) - FU  >*NU  +  2*U  *( - FX  >>#S 

dU  dX  2  2  1  2  dU  2  7 

1  2  dU  1 

1 

2  2 

d  d  d  d 

+  (-  U  *< - FX  >*S  +  < - FX  >*NU*S  +  ( - FX  -  2*( - 

1  dU  2  13  dU  2  14  2  2  dU  dX 

2  2  dX  12 

2 

d  d  d  d 

>  U  *< -  FX  )  -  U  *( -  FX  )  +  FU  J*S  +  U  *< -  FU  )*S  -  <-— 

2  dX  2  1  dX  2  2  7  1  dU  1  13  dU 

2  1  2  ; 


d  d  d 

-  ( -  FU  )*NU  +  U  *< -  FU  )  +  U  *< -  FU  >)*EPSLON 

2  1  2  dX  1  1  dX  1 

dX  2  1 


(088) 

<C89> 


BATCH  DONE 
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d 

*< - FX  > 

2  dU  2 
2 


FU  >>*NU 
1 


FU  >*NU*S 
1  14 


The  functions  can  be  simplified  considerably  by 

making  use  of  the  following  procedure:  if  all  the  independent 

coefficients  of  the  S  '  s  that  contain  only  a  single  term 

m 

are  set  to  zero,  these  conditions  on  the  groups  can  be  used 
immediately  to  discard  many  of  the  remaining  terms  of  5  . 

For  instance,  in  Hj  we  note  that 


FX, 


(A. 


’  fx2  *  %  fx2  ’  ijj  fx2  "  m-2  FOi 


By  using  these  relationships  many  of  the  remaining  terms  of 

can  be  crossed  out  immediately  from  the  printout.  For 
example,  by  using  Eqs.  (A.l)  in  Hx  the  terms  involving 
2 

S4  ,  S13S7,  and  S13  vanish.  By  equating  to  zero  the 
coefficients  of  the  remaining  terms  of  ,  those  involving 

S4,  S^,  and  the  constant  term,  we  obtain  the  following 
determining  equations 

aq  FX2 '  luj  ro2  -  iq;  FXi  +  roj  F0i  “  0  - 

-  hq  to2  +  Sbj-  fu2  *  0  - 
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